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Abstract 

We theoretically studied the in-plane tunneling spectroscopy of the hybrid structure composed 
of a metal and two-dimensional electron gas with Rashba spin-orbit coupling. We found that the 
energy spacing between two distinct features in the conductance spectrum can be used to measure 
the Rashba energy. We also considered the effect that varying the probability of spin-conserving 
and spin-flip scattering at the interface has on the overall conductance. Surprisingly, an increase in 
interface scattering probability can actually result in increased conductance under certain condi- 
tions. Particularly, in the tunneling regime, an increase in spin-flip scattering probability enhances 
the conductance. It is also found that the interfacial scattering greatly affects the spin polarization 
of the conductance in metal, but hardly affects that in the Rashba system. 
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I. INTRODUCTION 



Structural inversion asymmetry of the confining electrostatic potential results in an in- 
trinsic spin-orbit coupling of electrons in a two-dimensional (2D) electron gas (EG), which 



can be described by the Rashba Hamiltonian 



where p is 2D momentum, m* is the electron effective mass, j is the direction perpendicular 
to the plane of motion, A is the spin-orbit coupling parameter, which can be tuned by 
applying an external gate voltage perpendicular to the 2D plane, and the components of a 
are the Pauli spin matrices. The spin-orbit interaction lifts the spin degeneracy and causes 
the original parabolic energy spectrum to split into two branches: E^ ± = |^ ± h\k, where 
k is the magnitude of the wave vector. The density of states of this system is the same as 
that of the 2D free electron gas for all energies higher than the crossing point of the two 
branches. However, at the bottom of the band, the density of states has van Hove 
singularity because the minus branch has an annular minimum for k = k = m*X/h instead 
of a single-point minimum as in the free electron gas. These properties lead to interesting 



phenomena, like the spin 
spintronics (see e.g. Refs. 



rail effect (see e.g. Refs. |4j for review), and to applications in 
5I for review). 

The Rashba effect has been seen in many systems like surface alloys and semiconductor 
heterostructures. Several techniques have been used to study the spin-split states in these 
systems. For instance, angle-resolved photoemission spectroscopy , [hJ] and scan- 

ning tunneling microscopy llj are used in surface alloys. The former technique is utilized 
mainly to obtain the energy dispersion and the Fermi surface map, from which the effective 
mass, the magnitude of the band splitting, and hence the Rashba spin-orbit coupling en- 



ergy, E\ = S 2 ^/(2m*), can be extracted p, Q, Q . In the latter technique, the electric 
current is driven through a sharp tip perpendicular to the 2D plane and the differential 
conductance (dl/dV) spectrum can be obtained. One can deduce the Rashba energy by 
fitting the dl/dV spectrum to the local density of states of the 2DEG Hj. In both cases, to 
obtain information about the Rashba spin-orbit coupling, extensive data fitting is needed. 

In semiconductor heterostructures, the Rashba energy is measured using the Shubnikov- 
de Haas oscillations Q, [3]. The existence of the spin splitting at the Fermi energy leads to 
beating in the oscillations and the Rashba energy can be deduced from the position of the 



beating node. However, because this technique is done in the presence of magnetic field, it 
tends to provide an overestimate of the Rashba energy because it includes the effect of the 
Zeeman spin splitting jl^ ]. 

In this article, we propose a way to measure the spin-splitting energy more directly from 
the experimental data: the in-plane tunneling spectroscopy. In this technique, the Rashba 
energy can be determined by the energy difference between two features in the conductance 
spectrum. The required condition for the measurement is that the energy resolution of the 
tunneling spectra is at least of the order of the Rashba energy itself. This condition can be 
easily achieved in modern tunneling measurements [ijj]. 

An intriguing property of 2DEG with Rashba spin-orbit interaction is spin- dependent 
transport. Many theoretical investigations have shown that both electric and spin trans- 
port in hybrid structures between the Rashba system and various materials, like met- 



als 



18| , ferromagnets [18|, Il9j. 



20 



strength of the spin-orbit coupling [16 



2lfl . and superconductors [22j, are affected by the 



masses 
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20 



13, 



19 



20 



21 



221 ] , the inequality of the effective 



21| , and the transparency of the interface 19j, |20j, |22| . However, in these 



previous studies, only spin-conserving interfacial scattering was considered. In the pres- 
ence of interfacial spin-flip scattering, the equations describing the spin-up and spin-down 
spin states are coupled and one expects interesting consequences of this. For instance, in 
the study of the tunneling conductance spectrum of a semiconductor/superconductor junc- 
tion [231 ] . the non-spin-flip scattering, when present alone, is found to suppress the Andreev 
reflection process and hence the subgap conductance as expected. However, when the spin- 
flip potential scattering is also present at the interface, their combined effect surprisingly 
enhances the subgap conductance [23]. 

Here, we consider how the scattering potential barrier affects both the conductance spec- 
trum and the spin polarization of the conductance of a junction consisting of a metal and a 
Rashba system. As in previous work by Zutic and Das Sarma 23J], we find that the conduc- 
tance spectrum, which is usually suppressed in the presence of the interfacial scattering, can 
be enhanced by the combined effect of both types of scattering. We also find that the spin 
polarizations of conductance of the metal and the Rashba system are not equal. The spin 
polarization in the latter depends weakly on interfacial scattering, while that in the former is 
greatly affected. This suggests that a spin imbalance in the Rashba system is robust against 
variation in the quality of the junction interface. 
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FIG. 1: The top sketches are the energy contours of the electron in the metal (left) and the Rashba 
system (right). The angles 6 and tp are defined as those between the x axis and the momenta of 
electrons in the metal and the Rashba system respectively. The dashed line that crosses both sides 
shows the momentum states with the same k z . The dotted line is the line of the maximum value 
of k z , which defines the maximum incident angle m . The lower sketches are the corresponding 
energy spectra (E vs the magnitude of momentum). Ef and Eq are the metal Fermi energy and 
the off-set energy of the Rashba system respectively. 

In the next section, we describe the theoretical method and assumptions. In Section IHH 
we provide the results and discussion. The conclusion is given in the last section. 



II. METHOD OF CALCULATION AND ASSUMPTIONS 



We represent our junction by an infinite 2D system which lies on xz plane, where the metal 
and the Rashba system occupy the x < and x > region respectively. The two regions are 
separated by a flat interface at x — 0. The interfacial scattering is modeled by a Dirac delta 
function potential 24J. We consider ballistic transport in our junction. In the one-band 
effective-mass approximation, we describe our system by the following Hamiltonian: 

H = (pt-^P + V(x, z)) 1 + H R {x). (2) 



2m(x 

Each term is the 2x2 matrix acting on spinor states, p = —ih {p-jk + -2^)- The effective 
mass m(x) is position-dependent, i. e., [m(x)] _1 = m _1 @(— x) + (m*)~ l Q(x), where m and 
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m* are effective electron masses in the metal and the Rashba system respectively, and Q(x) 
is the Heaviside step function. V(x, z) is also a position-dependent function and is modeled 
by the expression 

V(x,z) = H5(x) + E Q(x) - E F Q(-x) (3) 

where H represents the scattering potential at the interface, E is the energy difference be- 
tween the Fermi level and the bottom of the plus branch (see FIG. [IJ, and Ep = h 2 q 2 F /{2m) 
is the Fermi energy of the metal. We assume that Ep is much larger than E . The diagonal 
elements of H, and correspond to the non-spin-flip scattering potential character- 
izing the quality of the junction, while = describe spin-flip scattering 23j. The 



Rashba Hamiltonian is written as 

H R (x) = 



25| 



3 



■ [\{x) (pxa) + (px a) \{x)\ 



(4) 



where X(x) = AG (a;). 

From the Hamiltonian, one can obtain the eigenstates and eigenenergy for the electrons 
in each region as follows. In the x < region, the energy spectrum is 

H 2 q 2 



E(q) 



2m 



E f 



(5) 



where q = a/ q 2 + q 2 is the magnitude of the 2D momentum of the electrons. In the x > 
region, the eigenenergy is obtained as 



k2 k 2 fc2r,2 



(6) 

2m* 2m* 

where k = >/k 2 + k 2 is the magnitude of the 2D momentum and ko = m*X/h. FIG. [1] shows 
the energy spectra and energy contours of the excitations in both sides of the junction. 

The wave function of the electrons in metal is written as a linear combination of incoming 
momentum state with equal spin components and a reflected state of the same energy and 




+ 



e -iq x x 1 e ik z z 



(7) 



where the 6f, are the amplitudes for reflection of spin-up and spin-down electrons respec- 
tively. q x = qcos8 and k z = gsin#, where 9 is the angle between q and the x axis. The 
magnitude of the momentum, q, depends on energy as 



q 



2m , 

-j^i-E + Ep) 



(8) 



Similarly, in the Rashba system, the wave function is obtained as a linear combination of 
two outgoing eigenstates of the same energy and k z \ 





COS ^ 


('■• 






isin^i 



sin 


2 


^ik x x 




cos 


2 


) 





(9) 



where (p^ are the angles between k and the x axis. For E > E , c + and c_ are the 
transmission amplitudes of electrons to plus and minus branch respectively. When E < E , 
c + and c_ refer to the transmission amplitudes of electrons to states with smaller and larger 
k of the minus branch respectively. The upper and lower signs in the first term of Eq. (J9j) 
are for E < E and E > E respectively, k^ = k^ cos v?fc± and k z = k^ sin tpk± , where the 
magnitudes of the momenta, k^, depend on energy as 



k- 
k~ 



^+\ik 2 + ^f(E-E ) 



± kn 



k 2 n + ^(E-E ) 



h 2 



(10) 

(ii) 



Again, in Eq. ffTTl) the upper and lower signs are for E < E and E > E respectively. The 
relationship between the angles (fk± and 9 is 



k^ sin (fk± = q sin 9. 



(12) 



We can obtain the probability amplitudes b^b^c + and c_ from the following matching 
conditions that ensure probability conservation [25]. 



m* dx 



^ M {x = 0,z) 



x=0 



dx 



x=0 



^ RS (x = 0,z) = * , 



Tfl 

2q F Z - i—k a z ) * , 

m* 



(13) 
(14) 



where Z = mH j [Y? qp) . The diagonal elements of Z will henceforth be referred to as 
Z u = Ztf and Zd = Z^, while the off-diagonal element will be denoted by Zp = Z^ = Z^. 
In what follows the spin flip term Zp will be responsible for the enhancement of a feature 
at the branch-crossing point in the conductance spectrum. 

The particle current density along the x direction is obtained from 

.„ 1 
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J. i 



&(x)v x fy(x) + (v x ^(x))^(x) 



(15) 
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where ^(x) is the spinor wave function, and v x = dx/dt = i [H(x),x] /h. From the current 
density, the reflection and transmission probabilities can be obtained: 



R^ — 


N 2 


Ri = 


IJ, |2 


T+ = 


m 
m* 


T_ = 


m 
m* 



q x 



(16) 
(17) 
(18) 

* J (19) 

where R-\,Ri are the reflection probabilities to spin- up and spin-down states, and T + ,T_ 

are the corresponding transmission probabilities. Also, the upper and lower signs in T + are 

for E < E and E > E respectively. As mentioned earlier, the matching conditions ensure 

that R 1 + R i +T + + T^ = 1. 

Since the current is independent of x, we consider the current density in the metal for 

simplicity. It can be written as a function of applied voltage V as follows. 

f x (eV)= ev x (l-R^R l )[f(E q -eV)-f(E q )} (20) 

g x >0,qz 

where e is the electron charge, v x is the x component of the electron group velocity, and 
f(E) is Fermi distribution function. 

By changing the integration variable and setting temperature to zero for simplicity, one 
can obtain the expression for the electric current as 



jtieV) = T~^r~~ dE d9 cos* Jl - — (1 - it! T - R L ) (21) 



h 2ix J Q J_ em V E F 

where C 2 is the area of the metal and 8 m is the maximum angle of the incident electrons 
from the metal (see FIG.[T]): 9 m = sin -1 (k~(E)/q(E)). Thus, the differential conductance 
G{V) = djl/dV at zero temperature is 



G ^ = t^F u M cos "v 1 - si (1 - flt - Ri) (22) 

The finite temperature will smear the features in the conductance spectrum but will not 
change their positions (assuming that the strength of the Rashba spin-orbit coupling does 
not depend on temperature). 

The spin polarization of the conductance V is defined as 

V(E) = g fc >o^(^T-n;t,J 

7 




FIG. 2: On the left side is the spectrum in the case where the energy band of the Rashba system is 
occupied (Eq = —0.075Ep) and on the right is that where the band is unoccupied (Eq = O.Q5Ep). 
Z and Zp are set equal to zero, m/m* = 10 and ko = 0.05qp. The derivative of the conductance 
spectrum on the right (dG/dV) is shown in the inset. 



where n a v a is the number of electrons with spin o that cross a plane of interest per unit time. 
In terms of the reflection probabilites this spin polarization in the metal can be written as 



Vm(E) 



E 



q x >0,q z 



E a 



'to>o,«, v 1 - R l - R i) 

and the spin polarization in the Rashba system in terms of the transmission probabilities is 
written as 

E g ,>o, , ( T + cos Vk+ - T_ cos cp k - ) 



(24) 



V R s(E) 



(25) 



As expected, the spin polarization measures the difference in number of the carriers with 
spin-up and spin-down on each side. 



III. RESULTS AND DISCUSSION 



In this section, we discuss the effect of the interfacial scattering on the differential con- 
ductance spectra and the spin polarization of conductance on each side of the junction. In 
all plots, we set m/m* = 10 and k = 0. 05(7 ??, which corresponds to typical experimental 



values in metal/Rashba system junctions 



23 



26| . In FIG. |2l two conductance plots for two 



values of E are shown. Positives values of E means the energy spectrum of the Rashba 



system is unoccupied and the positive eV across the junction will cause the current to flow 
from the metal to the Rashba system. 

In the case shown in FIG. HT1 where the energy spectrum is occupied (E = —0.075Ep), 
the results are identical in shape to those in the unoccupied case (E = +0.05Ep), but the 
applied voltage eV across the junction has to be negative. There are two main features at 
the voltage corresponding to the bottom and the branch-crossing of the energy band. The 
distance between them depends on E\, which is the quantity of interest. The value of Eq is 
not important, i. e., changing E causes a rigid shift in energy, and will henceforth be set 
equal to zero. 

We do not consider the spin filtering interface. That is, we set the non-spin-flip scattering 
strength Z u — Z& — Z. It is well-known that the difference in Z u and Z^ will cause a spin- 
filtering effect. That is, a higher Z u will make the transport of the spin- up electrons less 
favorable and vice versa. This effect cannot be seen in the conductance spectrum and will 
not be considered in the spin polarization. 

A. Differential conductance spectra 

In all conductance plots, the conductance is in units of e 2 C? qp j {2nK) . The conductance 
spectra G with different Zp in different limits of Z are shown in FIG. El Junctions with 
metallic contacts are characterized by Z <C 1, whereas those in the tunneling limit are 
characterized by Z > 1. In general, the conductance is zero until the applied voltage 
reaches eV = —Ex, which is the bottom of the band of the Rashba system. The conductance 
increases suddenly with large initial slope that decreases steadily until a second feature: the 
kink occurring at eV = 0, which is the crossing point of the two branches of the band. After 
this point, the conductance increases approximately linearly. In the presence of Zp, there 
occurs a discontinuity in the conductance at eV = 0. The height of the jump depends on 
both Z and Zp. The energy difference between the onset and the discontinuity in slope of 
the conductance spectrum can be used to measure the magnitude of the Rashba energy E\. 

In addition to the influence on the discontinuity at eV = 0, the interfacial scattering af- 
fects the overall conductance spectrum as well. For metallic contacts, the spin-flip scattering 
suppresses the conductance as expected. However, in the intermediate and the tunneling 
limits, the results are rather surprising. As can be seen in FIG. EJ^b) when Z = 0.5, the 
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FIG. 4: Differential conductance G(eV) plotted as a function of the spin-flip barrier height Zp at 
a constant energy eV slightly below [upper panel, denoted by G < {Zp)] and slightly above [lower 
panel, denoted by G > (Zp))] the energy corresponding to the crossing of the Rashba-split bands. 
The arrows indicate the values of Z F , where the maximum differential conductances G > and G < 
occur, for Z = 1.0 (thick arrows) and 2.0 (dashed-dotted arrows). 

increase in Zp from zero to a small value (less than 0.5) does not affect the conductance 
much. Only when Zp is increased beyond 0.5, does the conductance get suppressed. When 
Z is high, e. g. Z — 2.0 as in FIG. El^c), the conductance spectrum can be enhanced by the 
increase in Zp up to a value Z F , after which the spectrum becomes suppressed. Z F is found 
to depend strongly on Z. 

One can see the effect on the conductance spectrum of spin-flip scattering more clearly 
by considering plots of the conductance G as a function of Zp for energies just below and 
just above 0. In FIG. |G<e G{-5) and G> = G(+S), where 5/E x = 0.8, are plotted as 
a function of Zp for different values of Z. For small Z, both G > and G < decrease with Zp 
as expected. However, this trend starts to change when Z is higher than 0.5. That is, both 
G > and G < increase with Zp and reach a maximum value at Z F (as indicated by arrows in 
FIG. HJ, after which they decreases with Zp. Notice that Z F is a little smaller for G < than 
for G > and is approximately equal to Z. It should be noted that a similar dependence of 
both G > and G < on Z can also be seen, if one plots G > and G < as a function of Z instead. 
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FIG. 5: The plots of the spin polarization of conductance in metal and Rashba system (RS) as a 
function of energy when Z is (a) and (b) 0.5. 

B. Spin polarization of conductance 

The plots of the spin polarizations in both metal and Rashba system as a function of 
energy are shown in FIG. [5j The spin polarizations of the two sides are very different. In 
Rashba system, it is always negative, whereas in the metal it is positive when the spin-flip 
scattering is not strong. This may be understood by considering the density of states of the 
Rashba system. 

The density of states of the minus branch is larger than that of the plus branch. As we 
can see from FIG. |6j because the spins of the transmitted states of the minus branch are 
mostly pointing down, it is not surprising that the spin polarization in the Rashba system 
is negative. As for the metal side, because more spin-down states are transmitted into the 
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FIG. 6: Density of states of each branch of the 2DEG with the Rashba spin-orbit coupling. The 
contour plots on the left and on the right are those in the case where E < Eq and E > Eq 
respectively. When E > Eq, the outer contour is that of - branch and the inner one is that of + 
branch. When E < Eq, both energy contours belong to the - branch. The arrows represent the 
spin direction of the states with positive v x . 

Rashba system, the spin polarization is positive. 

The interfacial scattering does not affect the spin polarization in the Rashba system as 
much as in the metal. The increase in either Z or Zp seems to slightly change the magnitude 
of the spin polarization. However, in metal the interfacial scattering potential affects the 
spin polarization a great deal. For a particular value of Z, the increase in Zp can cause the 
spin polarization in metal to change sign. 

IV. CONCLUSIONS 

According to the results from our simple model, one can directly use in-plane tunneling 
conductance spectrum to measure the Rashba energy of a system with the Rashba spin- 
orbit coupling. The energy difference between the onset and the discontinuity in slope of 
the conductance spectrum is equal to the Rashba energy. Both features are found to be 
robust against variation in the quality of the junction. 

Experimentally, to be able to measure the Rashba energy, the required energy resolution 
is at least of the order of the Rashba energy itself and the temperature is low enough in order 
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that both features are visible. The Rashba energies in semiconductor-based heterostructures 



such as InAs, InGaAs and GaN, are of order 1 meV 
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321 ]. whereas those of 

oys like Li/W(110), Pb/Ag(lll), and Bi/Ag(lll) can be as large as 200 meV n[ 
351 ] . These conditions can be readily met in modern tunneling measurements 15]. 
We also found that as the current is driven through the system, an imbalance of spin in 
both sides occurs. The spin polarization of the conductance in the metal is found to depend 
strongly on both types of the interfacial scattering and can disappear when the barrier is 
in the tunneling regimes. On the contrary, in the Rashba system the spin polarization is 
always present and only slightly affected by interfacial scattering. This finding suggests that 
the spin imbalance caused by current flow in the system with the Rashba spin-orbit coupling 



is robust against variation in the quality of the junction as well. 
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